Disregarding multimappers leads to biases in the functional assessment of NGS data.

BACKGROUND: Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. RESULTS: In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. CONCLUSION: Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance.

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Integrative and Conjugative Elements and Prophage DNA as Carriers of Resistance Genes in Erysipelothrix rhusiopathiae Strains from Domestic Geese in Poland.

Goose erysipelas is a serious problem in waterfowl breeding in Poland. However, knowledge of the characteristics of Erysipelothrix rhusiopathiae strains causing this disease is limited. In this study, the antimicrobial susceptibility and serotypes of four E. rhusiopathiae strains from domestic geese were determined, and their whole-genome sequences (WGSs) were analyzed to detect resistance genes, integrative and conjugative elements (ICEs), and prophage DNA. Sequence type and the presence of resistance genes and transposons were compared with 363 publicly available E. rhusiopathiae strains, as well as 13 strains of other Erysipelothrix species. Four strains tested represented serotypes 2 and 5 and the MLST groups ST 4, 32, 242, and 243. Their assembled circular genomes ranged from 1.8 to 1.9 kb with a GC content of 36-37%; a small plasmid was detected in strain 1023. Strains 1023 and 267 were multidrug-resistant. The resistance genes detected in the genome of strain 1023 were erm47, tetM, and lsaE-lnuB-ant(6)-Ia-spw cluster, while strain 267 contained the tetM and ermB genes. Mutations in the gyrA gene were detected in both strains. The tetM gene was embedded in a Tn916-like transposon, which in strain 1023, together with the other resistance genes, was located on a large integrative and conjugative-like element of 130 kb designated as ICEEr1023. A minor integrative element of 74 kb was identified in strain 1012 (ICEEr1012). This work contributes to knowledge about the characteristics of E. rhusiopathiae bacteria and, for the first time, reveals the occurrence of erm47 and ermB resistance genes in strains of this species. Phage infection appears to be responsible for the introduction of the ermB gene into the genome of strain 267, while ICEs most likely play a key role in the spread of the other resistance genes identified in E. rhusiopathiae.

Stable, fluorescent markers for tracking synthetic communities and assembly dynamics.

BACKGROUND: After two decades of extensive microbiome research, the current forefront of scientific exploration involves moving beyond description and classification to uncovering the intricate mechanisms underlying the coalescence of microbial communities. Deciphering microbiome assembly has been technically challenging due to their vast microbial diversity but establishing a synthetic community (SynCom) serves as a key strategy in unravelling this process. Achieving absolute quantification is crucial for establishing causality in assembly dynamics. However, existing approaches are primarily designed to differentiate a specific group of microorganisms within a particular SynCom. RESULTS: To address this issue, we have developed the differential fluorescent marking (DFM) strategy, employing three distinguishable fluorescent proteins in single and double combinations. Building on the mini-Tn7 transposon, DFM capitalises on enhanced stability and broad applicability across diverse Proteobacteria species. The various DFM constructions are built using the pTn7-SCOUT plasmid family, enabling modular assembly, and facilitating the interchangeability of expression and antibiotic cassettes in a single reaction. DFM has no detrimental effects on fitness or community assembly dynamics, and through the application of flow cytometry, we successfully differentiated, quantified, and tracked a diverse six-member SynCom under various complex conditions like root rhizosphere showing a different colonisation assembly dynamic between pea and barley roots. CONCLUSIONS: DFM represents a powerful resource that eliminates dependence on sequencing and/or culturing, thereby opening new avenues for studying microbiome assembly. Video Abstract.

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution.

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety.

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Genome-Scale Investigation of the Regulation of azoR Expression in Escherichia coli Using Computational Analysis and Transposon Mutagenesis.

Bacterial azoreductases are enzymes that catalyze the reduction of ingested or industrial azo dyes. Although azoreductase genes have been well identified and characterized, the regulation of their expression has not been systematically investigated. To determine how different factors affect the expression of azoR, we extracted and analyzed transcriptional data from the Gene Expression Omnibus (GEO) resource, then confirmed computational predictions by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Results showed that azoR expression was lower with higher glucose concentration, agitation speed, and incubation temperature, but higher at higher culture densities. Co-expression and clustering analysis indicated ten genes with similar expression patterns to azoR: melA, tpx, yhbW, yciK, fdnG, fpr, nfsA, nfsB, rutF, and chrR (yieF). In parallel, constructing a random transposon library in E. coli K-12 and screening 4320 of its colonies for altered methyl red (MR)-decolorizing activity identified another set of seven genes potentially involved in azoR regulation. Among these genes, arsC, relA, plsY, and trmM were confirmed as potential azoR regulators based on the phenotypic decolorization activity of their transposon mutants, and the expression of arsC and relA was confirmed, by qRT-PCR, to significantly increase in E. coli K-12 in response to different MR concentrations. Finally, the significant decrease in azoR transcription upon transposon insertion in arsC and relA (as compared to its expression in wild-type E. coli) suggests their probable involvement in azoR regulation. In conclusion, combining in silico analysis and random transposon mutagenesis suggested a set of potential regulators of azoR in E. coli.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs.

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record.

BACKGROUND: Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes. RESULTS: Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses. CONCLUSIONS: We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Transposable elements regulate thymus development and function.

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/ß. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Molecular Epidemiology of Carbapenem-Resistant Klebsiella aerogenes in Japan.

Information regarding Klebsiella aerogenes haboring carbapenemase in Japan is limited. A comprehensive nationwide survey was conducted from September 2014 to December 2022, and 67 non-duplicate strains of carbapenem-resistant K. aerogenes were isolated from 57 healthcare facilities in Japan. Through genetic testing and whole-genome sequencing, six strains were found to possess carbapenemases, including imipenemase (IMP)-1, IMP-6, New Delhi metallo-ß-lactamase (NDM)-1, and NDM-5. The strain harboring blaNDM-5 was the novel strain ST709, which belongs to the clonal complex of the predominant ST4 in China. The novel integron containing blaIMP-1 featured the oxacillinase-101 gene, which is a previously unreported structure, with an IncN4 plasmid type. However, integrons found in the strains possessing blaIMP-6, which were the most commonly identified, matched those reported domestically in Klebsiella pneumoniae, suggesting the prevalence of identical integrons. Transposons containing blaNDM are similar or identical to the transposon structure of K. aerogenes harboring blaNDM-5 previously reported in Japan, suggesting that the same type of transposon could have been transmitted to K. aerogenes in Japan. This investigation analyzed mobile genetic elements, such as integrons and transposons, to understand the spread of carbapenemases, highlighting the growing challenge of carbapenem-resistant Enterobacterales in Japan and underscoring the critical need for ongoing surveillance to control these pathogens.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers.

BACKGROUND: While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts. RESULTS: Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family. CONCLUSIONS: This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria.

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery.

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome.

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

The evolution of transposable elements in Brachypodium distachyon is governed by purifying selection, while neutral and adaptive processes play a minor role.

Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.

Pangenome analysis reveals transposon-driven genome evolution in cotton.

BACKGROUND: Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution. RESULTS: Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes. CONCLUSIONS: Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Identification and characterization of endogenous retroviruses upon SARS-CoV-2 infection.

Endogenous retroviruses (ERVs) derived from the long terminal repeat (LTR) family of transposons constitute a significant portion of the mammalian genome, with origins tracing back to ancient viral infections. Despite comprising approximately 8% of the human genome, the specific role of ERVs in the pathogenesis of COVID-19 remains unclear. In this study, we conducted a genome-wide identification of ERVs in human peripheral blood mononuclear cells (hPBMCs) and primary lung epithelial cells from monkeys and mice, both infected and uninfected with SARS-CoV-2. We identified 405, 283, and 206 significantly up-regulated transposable elements (TEs) in hPBMCs, monkeys, and mice, respectively. This included 254, 119, 68, and 28 ERVs found in hPBMCs from severe and mild COVID-19 patients, monkeys, and transgenic mice expressing the human ACE2 receptor (hACE2) and infected with SARS-CoV-2. Furthermore, analysis using the Genomic Regions Enrichment of Annotations Tool (GREAT) revealed certain parental genomic sequences of these up-regulated ERVs in COVID-19 patients may be involved in various biological processes, including histone modification and viral replication. Of particular interest, we identified 210 ERVs specifically up-regulated in the severe COVID-19 group. The genes associated with these differentially expressed ERVs were enriched in processes such as immune response activation and histone modification. HERV1_I-int: ERV1:LTR and LTR7Y: ERV1:LTR were highlighted as potential biomarkers for evaluating the severity of COVID-19. Additionally, validation of our findings using RT-qPCR in Bone Marrow-Derived Macrophages (BMDMs) from mice infected by HSV-1 and VSV provided further support to our results. This study offers insights into the expression patterns and potential roles of ERVs following viral infection, providing a valuable resource for future studies on ERVs and their interaction with SARS-CoV-2.

Blessing or curse: how the epigenetic resolution of host-transposable element conflicts shapes their evolutionary dynamics.

Transposable elements (TEs) are selfish genetic elements whose antagonistic interactions with hosts represent a common genetic conflict in eukaryotes. To resolve this conflict, hosts have widely adopted epigenetic silencing that deposits repressive marks at TEs. However, this mechanism is imperfect and fails to fully halt TE replication. Furthermore, TE epigenetic silencing can inadvertently spread repressive marks to adjacent functional sequences, a phenomenon considered a 'curse' of this conflict resolution. Here, we used forward simulations to explore how TE epigenetic silencing and its harmful side effects shape the evolutionary dynamics of TEs and their hosts. Our findings reveal that epigenetic silencing allows TEs and their hosts to stably coexist under a wide range of conditions, because the underlying molecular mechanisms give rise to copy-number dependency of the strength of TE silencing. Interestingly, contrary to intuitive expectations that TE epigenetic silencing should evolve to be as strong as possible, we found a selective benefit for modifier alleles that weaken TE silencing under biologically feasible conditions. These results reveal that the dual nature of TE epigenetic silencing, with both positive and negative effects, complicates its evolutionary trajectory and makes it challenging to determine whether TE epigenetic silencing is a 'blessing' or a 'curse'.

Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae.

BACKGROUND: Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown. RESULTS: We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions. CONCLUSIONS: This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.